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We study the spin-s quantum Heisenberg antiferromagnetic models in a magnon theory free of 
any unphysical magnon state. Because the unphysical magnon states are completely removed in the 
magnon Hamiltonians and during the approximation process, we derive spin-s (s > i) reliable Neel 
temperature Tjv and reasonable sublattice magnetization unified for T < Tn- 
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Quantum Heisenberg antiferromagnetic (QHAFM) 
models are well-accepted models for insulating anti- 
ferromagnets. The insulating cupratcs La2C'u04 and 
RBa2Cu^OQ {R is an rare earth element) and their light- 
doped samples are typical antiferromagnets of spin 1/2. 
101 Because of their proximity to the high temperature 
superconductivity in the doped cuprates, they have ab- 
sorbed very much attention. On the other hand, insu- 
lating magnetic materials of high spin, such as the well- 
known Manganese materials, recently become important. 
H They can be described by the QHAFM models. But, 
for the QHAFM models, there has not yet existed any 
reasonable spin-s (s > i) Neel temperature Tjy and sub- 
lattice magnetization unified for all temperature T < Tjy. 

On the other hand, the defining spin operators in the 
QHAFM models are neither Bose nor Fermi operators. 
Many authors tried to transform the spin operators into 
the standard Bose (magnon) operators since Bloch in- 
vented the concept of spin wave, On the experimental 
side authors tend to make use of spin wave (magnon) 
theory to explain their experimental results. ||^ There- 
fore, magnon (or spin wave) theory is a very powerful 
approach to the QHAFM models. In this approach one 
makes use of Holstcin-Primakoff Q or Dyson trans- 
formation to transform the original spin Hamiltonians 
into magnon Hamiltonians. On a single site, the magnon 
state space is infinitely dimensional, but there are only 
2s -|- 1 physical spin states. There exist infinite extra un- 
physical magnon states in any magnon Hamiltonian. In 
removing effect of the unphysical states, many methods 
were proposed, J/j-p^ but the problem of the unphysical 
magnon states remained unsolved. In a magnon theory 
free of the unphysical states, ||l^] all of the unphysical 
states were put infinitely high in energy and their effect 
was removed in the magnon Hamiltonians and during the 
approximation process. In the case of spin 1/2, reliable 
ferromagnetic magnetization and antiferromagnetic sub- 
lattice magnetization were obtained under a simple but 
sensible approximation. 

In this Letter we study the (s > 1) quantum Heisen- 
berg antiferromagnetic models in the magnon theory free 
of the unphysical states. In the way similar to those in the 
previous works, ||l3|| we remove effect of the unphysical 
magnon states in the magnon Hamiltonians and during 



the approximation process, and thereby derive for the 
first time reliable analytical Neel temperature Tn and 
sublattice magnetization, being reasonable from zero to 
T/v, of the quantum Heisenberg antiferromagnetic model 
of any spin (s > 1) in three dimensions. The zero- 
temperature normalized sublattice magnetization Mq/s 
increases monotonously with spin s. The reduced sublat- 
tice magnetization M/Mq decreases monotonously with 
spin s increasing at any given reduced temperature T/Tn 
between and 1. Our sublattice magnetizations are sub- 
stantial improvements on those of the existing nonlinear 
spin wave (magnon) theories, which became unreason- 
able at high temperature. Our Neel temperature in the 
spin-i case is only 4% larger than the series expansion 
resuh. 

The standard spin-s QHAFM Hamiltonian is defined 

by 



-s+) + St 
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The summation is over the nearest neighbor sites of the 
lattice. The exchange constant Jy is positive and, in 
general, is direction dependent to describe real space 
anisotropy. We suppose our lattice can be divided into 
two sublattices so that the nearest neighbors of any site 
of one sublattice belong to the other sublattice. Making a 
TT rotation in the spin space of one of the two sublattices, 
the Hamiltonian (1) becomes 



(2) 



To transform the spin operators into the magnon opera- 
tors. We choose Holstein-Primakoff (HP) transformation: 



= aj\/2s - Hi, 5*+ = y/2s - rnai, SI ^ s - n^, (3) 

where the magnon number operator Ui is defined by 
rii = a\ai. The magnon operators a are standard Bose 
operators. After substituting the magnon transformation 
(y) into the Hamiltonian (||), we obtain the HP magnon 
Hamiltonian: 
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TT this reason, we construct the two series of Green func- 

fm fm — 1 



(4) 



where Hq — —^JZs^N , e = JZs, and Ai = -/2s — n^. 

The magnon Hamiltonian is different from the orig- 
inal spin Hamiltonian because the magnon operator 
on a single site has infinite states. For spin s we have 
only 2s + f physical spin states on a single site. To re- 
move effect of the infinite unphysical magnon states, we 
introduce an infinite-[/ term: 
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into the HP magnon Hamiltonian. As a result, our total 
Hamiltonian is defined by H' = Hm + Hu- The extra un- 
physical states are raised infinitely high in energy by the 
U term and actually decouple from the 2s + f physical 
states. The magnon square root Ai = y/2s — rii in the 
Hamiltonian can be expanded in terms of a|™a™. The U 
term makes the unphysical onsite correlation ^a|™a 

(m > 2s + 1) exactly equivalent to zero. So that the ex- 
pansion of \/2s — Hi is accurately truncated into a sum 
of 2s -|- 1 magnon operator product terms: 
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Substituting the Ai expression (|^) into the magnon 
Hamiltonian (^), our total Hamiltonian becomes a finite 
sum of rational magnon operator products. Therefore, 
the infinite-C/ term Hu removes effect of the unphysical 
magnon states and at the same time simplify the magnon 
Hamiltonian. |Q 

Being similar to the previous work, we work in the 
framework of Green function method and the approach of 
equation of motion. We define 1^{X\Y)) to be the Green 
function of the operators X and Y in Zubarev notation. 
The equation of motion of operator X is defined to be 
ij-^X = \X,H% or zX = [X,H'] in the frequency (z) 
space. Since the magnon number is not conserved for 
the AFM models, we need two kinds of Green functions. 
Intuitively, we have ((ai|a^)) and ((a-|a])). But they 
cannot compose a closed set of equations of motion be- 
cause of the infinite U term in the magnon Hamiltonian. 
On-site multimagnon Green functions inevitably appear 
due to the U term. We need a series of the on-site mul- 
timagnon Green functions to construct some composite 
Green functions so that their equations of motion are 
closed and contain terms of at most 1/U order. For 



tions: {{B,\Dp) and {{Bj\Dp) (to 



where 



and Bi ~ Ai 



1,2,3,- • -,23), 
Here the opera- 



tor Ai is equivalent to Ai with the al'^^af'' term removed. 
The Green functions satisfy the equations of motion 

(z - e)m\Df)) = WrS,, + J2 M{Tu\Df)) (8) 



and 

(z + e){{Bl\D]^)) = W^k, - M{Tl\Df)), (9) 

where = {[A,a„D^]), = l\a\A\,D';^\) , and 

Til = {s — ni)alAi — niAiUi. In Eqs. (||) and (^), other 
terms of 1/J7 order or smaller terms are neglected since 
they are zero when U tends to infinity. We make the 
on-site approximation riiXi = {ni)Xi ~ {n)Xi for the 
operator riiXi only ii i ^ I. Under this approximation, 
the equations of motion reduce to 

[z - ^)m\Df)) ^ wrk, + E '^^'(^ - ^)((^/i^r)) 

(10) 



and 



(z 7)((Sl|i?f )) =-Y^J.i{s- n)m\DJ^)), (11) 



where 7 = JZ^s — n). The parameter J is the largest one 
of the exchange Jij and the effective coordination number 
Z is defined hy Z = J2j(i) Jij/Jj where the j summation 
is made over all nearest neighbor sites of site i. It is 
not difficult to solve the above equations by transform- 
ing the Green functions into the fc-space. The resultant 
Green functions in the fc-space read: 



and 



W"^[z + JZ{s-n)] 
z2 - J2z2(s-n)^{l-rl) 
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(12) 



(13) 



The spectrum function is defined by = 
{"2/ JZ)J2i(j) Jjiexp{iRi — iRj). It is obvious that the 
magnon correlation functions {D™Bi) = V^^ and WJ" 
both are site independent. Therefore we obtain the self- 
consistent equations 



and 



^ iV ^ 2 



- 1 



(14) 



(15) 



2 



For spin s, we have only 2s non-zero on-site correlation 
functions: p„ = (Pni) ^ n < 2s), where the on-site 
multimagnon operator Pni is defined by Pni — ^a'^a'^. 
The unphysical correlation functions P„ (n>2s-\- 1) are 
exactly equivalent to zero. The self-consistent equation 
( p^ ) can be solved by 



2S-I-1 



7p' 



2s+l 



2s— n 

k=0 



(16) 



where n < 2s and the coefficient C," = nl/[ml{n — m)!] 
is the binomial coefficient. Since — s ~ rii and n = pi, 
we derive the following spin magnetization: 



{S')^s-ij + 



i2s + l)ijj 



2s+l 



2s+l 



2s+l 



(17) 



At very low temperature, ip can be expanded as ip = 
ipo + cT"^ (c is a positive constant . ) , where t^o = (C ~ 1 ) /2 
and ^ = l/-\/l — r^. As a result, we derive 

(S"^) = S'q — ciT^ (ci is also a positive constant.). Sq 



is given by the right-hand side of formula ([1^) with ip 
replaced by ipQ. The Neel temperature Tjy is given by 

JZ 

Tn^—s{s + 1), (18) 
3?7 

where the constant 7y is given by rj = X^fc 1/(1 ~ ^fc)- 
When temperature approaches to the Neel temperature 
Tjv, the sublattice magnetization {S^) is given by the 
asymptotic expression: 



(19) 



The positive constant a is defined by a = |s(s + 1)[1 + 
I -I- ]^(2s^ -|- 2s — 9)]~^/^. In one dimension, the integra- 
tions r/ and ^ both diverge. There is no antiferromagnetic 
order at any temperature. In two dimensions the integra- 
tion ^ converges, but the integration ry diverges. There is 
antiferromagnetic order only at zero temperature. These 
are consistent with the Mermin- Wagner theorem. Jisj 

At general temperature, one cannot derive any analyt- 
ical expression. We present some digital results in three- 
dimensional (3D) simple-cubic lattice. In Fig. 1 we show 
the sublattice magnetizations of spins 1/2, 1, 3/2, and 2 
as functions of temperature. For convenience, we normal- 
ize the saturation magnetizations into one, that is, we di- 
vide the magnetizations by the corresponding spin s. The 
normalized sublattice magnetizations increase with spin 
s increasing at zero temperature. For comparison with 
the results of the existing spin wave (magnon) theories, 
we also present the sublattice magnetizations of a typi- 
cal version of the existing nonlinear spin wave (magnon) 
theories (NLSW). It is obvious that the NLSW magneti- 
zations are unreasonable at high temperature, especially 
near the Neel temperatures. In fact, all of the exist- 
ing NLSW theories yielded double- valued magnetizations 



at high temperature. ||l^,|r^ Therefore the existing spin 
wave (magnon) theories work at most at low tempera- 
ture, and our magnon theory free of unphysical states 
works well in whole temperature region. 
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FIG. 1. Sublattice magnetizations (divided by s) as func- 
tions of temperature (in unit of J) for spin 1/2 (see Ref 13), 
1, 3/2, and 2. The solid lines are the results of the magnon 
theory free of any unphysical states. The dashed lines are 
results from other magnon (spin wave) theories (Ref 16) . The 
latter are not reasonable at high temperature (Refs 16, 17). 
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FIG. 2. The reduced sublattice magnetizations (M/Mq) as 
functions of reduced temperature {T/Tm) for spin 1/2 (see 
Ref 13), 1, 2, 4, 9, and oo. The reduced sublattice magne- 
tization decreases with reduced temperature increasing for a 
given spin s. It also decreases with spin s increasing at a given 
reduced temperature. 

To compare the magnetizations of different spins, we 
show the reduced sublattice magnetizations M/Mq of 
spins 1/2, 1, 2, 4, 9, and oo as functions of the re- 
duced temperature T/T^. Near the Neel temperature, 
the reduced magnetization behaves as \/l — T/Tn for 
any spin s. At low temperature, the reduced sublat- 
tice magnetization can be expanded as 1 — C2T^ (c2 is 
a positive constant.). For very large spin s, the be- 
havior holds within a temperature region of the width 
s. But the Neel temperature is proportional to s^. As 
a function of the reduced temperature t = T/Tj^, the 
the t region of square behavior shrinks to zero when spin 
s tends to infinity. Therefore, for spin oo, the reduced 
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magnetization is linear in the reduced temperature at 
T/Tm = t ^ Q. It is very interesting that the reduced 
magnetization monotonously decreases with spin s at any 
given reduced temperature t between and 1. 

The above results are obtained under the postulation 
of isotropy ( = J) in the 3D simple cubic lattice. In 
fact, the approach and the on-site approximation can be 
successfully applied to the cases of anisotropic exchange 
constants where Jy is direction dependent. In the cases 
of the cuprates, one has an in-plane exchange constant 
J and an interplane exchange constant 5 J {5 is small.). 
|p],p^ For line-like materials, one has an on-line exchange 
constant J and an interline constant 5J with small 5. 

The method and approximation process we make use 
of in this paper are similar to the previous work on the 
quantum Heisenberg spin-i ferromagnetic models and 
spin-i antiferromagnetic model. |Q In those cases, our 
magnetizations and sublattice magnetization were equiv- 
alent to the reasonable results which had been obtained 
by means of the direct spin-operator methods. [l9|-]2l|] In 
the case of spin 1/2, we obtain Tn = 0.989 J in the simple 
cubic lattice, being only 4% larger than the series expan- 
sion result Tn = 0.951J. Therefore, the rehabihty 
of our spin-s sublattice magnetizations and Neel temper- 
atures in this paper can be justified by the equivalence 
in the ferromagnetic cases and the comparison with the 
accurate series expansion result, jlj] 

In summary, we have studied the quantum Heisenberg 
antiferromagnetic models of spin s in a magnon theory 
free of the unphysical states. Effect of the unphysical 
states is completely removed in the magnon Hamiltonians 
and during the approximation process. We thereby de- 
rive reliable analytical the Nccl temperatures and sublat- 
tice magnetizations unified for all temperatures T < 
and obtain the multimagnon correlation functions for any 
spin. At any given reduced temperature, the reduced 
sublattice magnetization monotonously decreases with 
spin s increasing. Our sublattice magnetizations are sub- 
stantial improvements on those of the existing spin wave 
(magnon) theories. Our Neel temperature in the case 
of spin 1/2 is only 4% larger than the series expansion 
result. 
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